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We study the breakdown of the Stokes-Einstein (SE) and Debye-Stokes-Einstein (DSE) relations 
for translational and rotational motion in a prototypical model of a network-forming liquid, the 
ST2 model of water. We find that the emergence of "fractional" SE and DSE relations at low 
temperature is ubiquitous in this system, with exponents that vary little over a range of distinct 
physical regimes. We also show that the same fractional SE relation is obeyed by both mobile and 
immobile dynamical heterogeneities of the liquid. 

PACS numbers: 



The translational (D) and rotational (D r ) diffusion 
constants of a macroscopic object in a simple liquid are 
well-characterized by the Stokes-Einstein (SE) relation, 
Dt/T = c, and the Debye-Stokes-Einstein (DSE) rela- 
tion, D r r/T = c r . In these relations, r is a relaxation 
time proportional to the viscosity of the liquid, T is the 
temperature, and the value of the constants c and c r 
depend on the geometry of the object and the bound- 
ary conditions. While these relations were originally for- 
mulated for the diffusion of macroscopic objects, they 
are known to hold when the diffusing object is itself a 
molecule of the liquid 0, 0- Hence the SE and DSE 
relations provide a simple connection between mass and 
momentum transport in a liquid. 

As T decreases and the liquid becomes increasingly 
viscous (assuming crystallization is avoided), numerous 
experiments have shown a failure of the SE relation for 
T < 1-3TV, where T g is the glass transition tempera- 
ture dH IE IS 0,11 Specifically the ratio Dt/T is found 
to increase by as much as 2 or 3 orders of magnitude on 
cooling 8]. Translational diffusion is thus said to be en- 
hanced relative to viscosity. When the SE relation fails, 
it has been empirically found that a "fractional" Stokes- 
Einstein (F-SE) relation, D ~(t/T) _ £ holds for a wide 
range of liquids 0, H ES , where 0.5 < £ < 0.95, 
with values between 0.7 and 0.8 being most commonly 
reported. While a definitive theoretical prediction of the 
value of £ has been elusive, various theoretical models 
predict £ in the same range 0, 0, Q, . The situa- 
tion is more complex for the DSE equation. In the same 
T range where the SE equation fails, the DSE equation 
has been experimentally found to be valid for most liq- 
uids However, in other experiments H, 0, a 
fractional Debye-Stokes-Einstein (F-DSE) relation, D r ~ 
(r/T)~* r , has been observed. Further complicating the 
situation is the fact that D r is experimentally inacces- 
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sible, and experiments must test the DSE relation us- 
ing a re-orientational correlation time, which may not 
be simply related to D r . Molecular dynamics simula- 
tions can directly evaluate D r , and simulations of dumb- 
bell molecules show a failure of the DSE relation using 

A commonly proposed explanation for the breakdown 
of the SE relation is the presence of dynamical hetero- 
geneity (DH) 0,113, i.e. spatially correlated regions of 
relatively high or low mobility that persist for a finite life- 
time in the liquid, and that grow in size as T decreases. 
The existence of DH has been confirmed and quantified 
in both experiments 0, and simulations [2J. How- 
ever, only a few liquid simulation studies directly probe 
the relationship between DH and the breakdown of the 
SE relation [HEl. 

In this Letter, we study the breakdown of the SE and 
DSE relations in simulations of the ST2 model of wa- 
ter j24j . The ST2 model is a tetrahedral arrangement of 
charges, centered within a Lennard-Jones envelope, that 
qualitatively reproduces a number of water's anomalies. 
However, our goal is not to elucidate the properties of wa- 
ter specifically; indeed there are many water models more 
accurate than ST2. Rather, we choose this model be- 
cause it is known to display a diversity of extreme liquid 
behavior within a single system (Fig.^l. At high density 
p, ST2 behaves as a typical fragile glass-forming molecu- 
lar liquid. However, near p = 0.93 g/cm 3 and T — 240 K, 
a liquid-liquid critical point occurs in the equation of 
state, demarcating the onset of phase separation between 
a low density liquid (LDL) and a high density liquid 
(HDL) phase |2g. At still lower p, near 0.83 gm/cm 3 , 
the behavior of ST2 is dominated by the emergence of an 
increasingly well-structured random tetrahedral network 
(RTN) as T decreases. The emergence of the RTN is 
associated with the onset of a fragile-to-strong crossover 
in the liquid transport properties |2f| . We study each of 
these three distinct regimes, to test how the character of 
SE and DSE breakdown changes as the behavior of the 
liquid state itself changes. We also test if the origins of 
the F-SE relation can be found in the DH of the system. 

Unless otherwise indicated, our data are generated 
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FIG. 1: ST2 phase behavior in the p-T plane. Solid curves 
show spinodals bounding regions of liquid-gas and HDL-LDL 
liquid phase separation. C is the critical point of the HDL- 
LDL transition. The three arrows indicate the densities fo- 
cused on in subsequent plots. The dotted line is the locus of 
density extrema. The dashed curve is the locus of points at 
which Dt/T is 1.5 times its value at 400 K, along isochores. 



by molecular dynamic simulations using TV = 1728 
molecules with a time step of 1 fs. Simulations consist 
of equilibration, followed by a production phase, each of 
which is run for the longer of 100 ps or the time needed 
for the mean-squared displacement per molecule (r 2 ) to 
reach 1 nm 2 - roughly three molecular diameters. Simu- 
lations are carried out in the range 250 K< T < 400 K at 
intervals of 5 K and for 0.80 g/cm 3 < p < 1.20 g/cm 3 at 
intervals of 0.01 g/cm 3 . For simplicity of the figures, we 
only present data along three representative isochores: 

(i) p = 0.83 g/cm 3 , where the RTN structure dominates; 

(ii) p — 0.93 g/cm 3 , approximately the critical density for 
the liquid- liquid transition; and (hi) p — 1.03 g/cm 3 at 
which the system behaves more as a simple fragile liquid. 
The simulations are carried out at fixed p and T, employ- 
ing the Berendsen heat bath with a time constant of 2 ps. 
Electrostatic interactions are truncated at 0.78 nm and 
the energy and pressure are corrected using the reaction 
field method %2$. 

We evaluate D from the asymptotic behavior of (r 2 ) = 
6Dt; and D r by tracking the mean-squared angular dis- 
placement (4> 2 ) of the molecular dipole moment vector 
and use the asymptotic relation (0 2 ) = 4D r t j2S|. The 
alpha-relaxation time r is defined as the time where the 
normalized coherent intermediate-scattering function de- 
cays to a value of 1/e, evaluated at the closest discrete 
wave number to q = 18 nm , the approximate loca- 
tion of the first peak in the static structure factor. The 
dipole relaxation time t-i is defined as the time at which 
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FIG. 2: (a) Breakdown of the SE relation at low T. (b) Test 
of the DSE relation using D r . (c) Test of the DSE relation 
using 1/t2 as a proxy for D r . 



the correlation function (P2[cos9(t)]) decays to a value 
of 1/e, where P2 is the second Lengendre polynomial and 
9{t) is the angle defined by the orientations of the dipole 
moment vector of a molecule at times t and zero. 

To determine in which regions of the phase diagram 
the SE and DSE relations are valid, we show the be- 
havior of Dt/T and D r r/T in Fig. For T > 300 K, 
Figs. Efa) and (b) show that both the SE and DSE ra- 
tios are nearly constant, as expected. The growth of the 
SE relation reflects the expected enhanced translational 
diffusion. In the same sense, we also find enhanced rota- 
tional diffusion; noting that Fig.^b) employs a log scale, 
the enhancement of rotational motion exceeds that of 
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FIG. 3: Testing for (a) F-SE behavior; (b) F-DSE behavior 
using D r ; and (c) F-DSE behavior using \jr%. Curves for 
p — 0.93 and 1.03 g/cm 3 have been multiplied by arbitrary 
factors to facilitate comparison. Solid lines have the slopes 
indicated. 



translational motion. This is consistent with spin-lattice 
relaxation experiments on water j^. The T at which 
the SE relation breaks down is weakly p dependent, as 
shown in Fig. ^ 

The deviation from the DSE relation can be seen as 
somewhat surprising since many (but not all) experi- 
ments have found that the DSE relation appears to be 
valid even when the SE relationship fails 0, 0| . In experi- 
mental studies, D r cannot be directly measured. As a re- 
sult, the inverse dipole relaxation time t^ 1 is frequently 
used as a substitute, since it can be readily measured. 
Since we found a deviation from the DSE relation using 
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FIG. 4: Testing for F-SE behavior at p = 0.83 g/cm 3 . Circles 
are data from N = 1728 simulations. The "clouds" of other 
symbols are obtained from ensemble simulations of a N = 216 
system at 255 K (+), 260 K (y), 265 K (o), 270 K (A) and 
275 K (□). The "cloud" for each T has been repeated and 
shifted by an arbitrary factor to reveal its shape in isolation 
from the other clouds. Solid lines have the slopes indicated. 



D r , we consider whether the using T2 yields different re- 
sults, as shown in Fig. Hjc). Curiously, using we find 
no T range where the DSE relation appears valid. For 
the system we study, the difference found using T2 versus 
D r may result from large reorientational motions associ- 
ated with hydrogen bond switching. Whatever the cause, 
Fig. |2£c) illustrates that substitution of D r with T2 may 
not, in general, be appropriate [l6| . 

To test if the SE relation is replaced by the F-SE form 
at low T, we parametrically plot the relation between be- 
tween D and t/T in Fig.jSfa). We find that there are two 
distinct regions: (i) T > 275 where we find an exponent 
£ = 1, consistent with the SE relation, and (ii) lower T 
where we find a F-SE relation with an exponent in the 
commonly observed range 0.7 < £ < 0.8. Moreover, the 
value of £ is nearly independent of p. Since the DSE 
relation is also clearly violated in this system, we check 
for a F-DSE relation in Fig. EJb). Indeed, analogous to 
Fig. EH a), we find that a high-T DSE regime is replaced 
at low T by a F-DSE relation with £ r s=i 0.25, signifi- 
cantly smaller than £ for the F-SE relation. The smaller 
value of £ r compared to £ is also observed in experiments 
on OTP [16j and is consistent with our finding that the 
rotational motion is more dramatically enhanced than 
the translational diffusion. Like £, we find that £ r varies 
weakly with p, if at all. Hence, despite the wide variation 
of liquid state properties over the density range studied 
here, the quantitative form of the F-SE and F-DSE be- 
haviors is strikingly uniform. 

We next analyze the behavior of the ST2 model to test 
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for a connection between DH and F-SE behavior. We 
focus on the p = 0.83 g/cm 3 isochore, and conduct 40 
independent simulations with N — 216 molecules at each 
T for 255 < T < 275 K, at 5 K intervals, again using 
the criterion that each run continues until (r 2 ) = 1 nm 2 . 
The smaller system (compared to N = 1728) allows us to 
probe much longer time scales and to generate an ensem- 
ble of many runs in a computationally accessible time. As 
DH emerges at low T, each member of the ensemble at 
the same T becomes more strongly influenced by tran- 
sient mobile and immobile regions. As a consequence, 
some members of the ensemble take a longer time to reach 
the run-time criterion, and some shorter. Thus, we can 
study the differences between systems dominated by mo- 
bile and immobile DH subdomains without the need to 
classify individual molecules by their mobility. 

It has been hypothesized that the SE violation is 
due primarily to the presence of transient mobile re- 
gions 0, 0, |23. In this picture, a "background" of 
largely immobile regions obeys the normal SE behavior, 
while mobile regions violate the SE relation. If true, we 
should find that the ensembles members that are domi- 
nated by mobile regions should exhibit a prominent F-SE 
behavior, while the members dominated by immobile re- 
gions should tend to follow a normal SE behavior. In this 
case, a parametric scatter plot of D against t/T (Fig.0J| 
for all ensembles should show a local curvature of the 
data for each ensemble at a given T, such that the low-D 
envelope of the data reflects SE behavior, and the high- 
D envelope gives F-SE behavior. This is not what Fig.0] 
shows. Rather, we find that all systems, both those that 
are relatively fast and those that are slow, all obey the 
same F-SE relation with £ rj 0.77. Therefore, for the 
ST2 model of water at p = 0.83 g/cm 3 , both the mobile 
and immobile regions of the liquid deviate in the same 
way from the SE relation. In this sense, the violation 
of the SE relation is spatially homogeneous, despite the 
existence of DH in the system. 

We note that £ = 0.77 matches the prediction of the 
entropic barrier hopping theory of Ref. |14|. although 
Ref. [14( assumes the SE relation is valid over small do- 
mains, and that FSE behavior arises from averaging over 
spatial variation of the domains. However, more recent 



work by the same authors [30j observes a spatially ubiq- 
uitous FSE, as found in this work. Very recently, Ref. [22] 
studied the impact of heterogeneous dynamics on the SE 
relation in the hard sphere system by identifying a sub- 
set of "hopping" mobile molecules; their results indicated 
that only the hopping molecules violate the SE relation, 
though they did not investigate F-SE behavior. It would 
be useful to conduct a similar analysis of our system, in 
which hopping can also be expected at low T. 

Overall, our results demonstrate the remarkable ro- 
bustness of the F-SE and F-DSE behaviors to variations 
in both the mean and local molecular environment. The 
fractional behavior is observed across three distinct phys- 
ical regimes (fragile, critical, and network- forming), with 
little variation of the exponents £ and £ r , indicating an 
almost complete insensitivity to changes in the average 
liquid structure. Even the large static fluctuations asso- 
ciated with the approach to a second order critical point 
do not observably effect the manifestation of the F-SE 
and F-DSE relations. We also find no sensitivity of the 
F-SE behavior on the relative proportion of mobile and 
immobile dynamical heterogeneities within a given sys- 
tem. Both locally fast and slow regions of the system 
obey the same F-SE relation, so the phenomenon cannot 
be attributed to any one extreme subset. We emphasize 
that our results do not imply that F-SE behavior and 
DH arc disconnected, since the two phenomenon emerge 
in the same range of T 31], only that the relationship 
is not the anticipated one. Our analysis also does not 
indicate a specific alternative origin for F-SE and F-DSE 
behavior. However, it does suggest that in a success- 
ful theory of the phenomena, F-SE and F-DSE behavior 
should emerge as an intrinsic behavior of the liquid state 
as T — > T g , insensitive to the details of the structural and 
dynamical fluctuations occurring in the system. This is 
certainly consistent with the nearly ubiquitous observa- 
tion of SE breakdown in glass forming liquids. 
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